Can X(3872) be a J p = 2 tetraquark state? 
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In this article, we test the nature of V(3872), which is assumed to be a P-wave 
[cg]-scalar-diquark [cg]-axial-vector-antidiquark tetraquark state with J p = 2~. The 
interpolating current representing the J p = 2~ state is proposed. Technically, con- 
tributions of the operators up to dimension six are included in the operator product 
expansion (OPE). The mass obtained for such state is m 2 - = (4.38 ± 0.15) GeV. 
P-f We conclude that it is impossible to describe the X(3872) structure as J p = 2~ 



tetraquark state. 

PACS numbers: 11.55.Hx, 12.38.Lg, 12.39.Mk 
The state AT (3872) was first discovered by Belle 



l| in the ti + -k J ftp mode and then 
confirmed by the CDF [2], D0 [3], and BABAR \^\ Collaborations in the same decay 



< 
(N 

9-; 

channal. The most recent measure of its mass is [5[ 

<N ; m x = (3871.85 ± 0.27(stat) ± 0.19(syst)) MeV, (1) 

. with a width of T x < 1.2 MeV. Belle [lj and CDF [6| propose that it proceeds through 

the X — > J/tjjp — » J/ifm +/ K~ decay. Since a charmonium state has isospin zero, it can 
not decay into X — > J/ipp, so the X(3872) is identified as an "exotic" state. According 
to the CDF analysis of the decay angular distribution [6] and the invariant -k + tt~ mass 
distribution [7] of the J/ip7r + 7r~ decay mode, only 1 + and 2~ assignments are possible. 
The close proximity of X(3872) mass to the DD* threshold indicates that X(3872) might 
^ ' be a loosely bound DD* molecular state, whose quantum number is J p = 1 + . Also, an 

angular analysis applied to the 2tc mass distribution in J/ipp favors the quantum number 
J p = 1 + [sj]. In compliance with these quantum numbers, many literatures have appeared 
in the past years. Its possible interpretations include the molecular state, tetraquark state 
and hybrid charmonium (see reviews [^-[lll and references therein). Using QCD sum 
rules (QCDSR) [3], Nielsen et al. discuss the possibility that it is possible to describe the 
X(3872) structure as a mixed molecule-charmonium state and study its strong decay and 
radiative decay 

HQ. 

Very recently, the BABAR collaboration has performed angular distribution analysis of 
the decay B — > J/ipuK, indicating that P-wave between J/ip and u is favored, so that 



quantum numbers J p = 2~ is preferred 19|. In this case, the most conventional explana- 
tion is the l 1 /^ charmonium state rj C 2(lD). In Ref 20J], the radiative transition processes 
Tjc2(lD) — > J/tp(ip') + 7 is investigated within several phenomenological potential models 
with the assumption that V(3872) is a r] c2 (l-D) charmonium, which are in contradiction 
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with the existing BABAR measurements 21] . The data on its D D°tt decay mode |22j 



also contradict the l 1 /^ charmonium interpretation of the A(3872) 23]. The decay of 



B — > rj&X is studied in NRQCD factorization framework, which indicates that X(3872) 



is unlikely to be a l l D 2 charmonium state 24|. Thus, we have to resort to exotic expla- 



nations for the J p = 2" quantum numbers. In Ref. 25], it's shown that the molecular 
interpretation appears to be untenable, but the tetraquark interpretation may be a vi- 
able candidates to be A (3872) with J p = 2~ . Follow their opinion, we study the mass 
of A (3872) as a P-wave [cg]-scalar-diquark [eg] -axial- vector- antidiquark tetraquark state 
with J p = 2- using the QCDSR. 

The interpolating current representing a J p = 2~ P-wave tetraquark state with [cq\- 
scalar-diquark and [eg] -axial- vector- antidiquark fields is adopted as 

= "-^[{qlC^D^Ccl )] . (2) 

Herein the index T represents matrix transposition, C means the charge conjugation ma- 
trix, denotes the covariant derivative, while a, b, c, d, and e are color indices. 

In the QCDSR approach, the mass of the particle can be determined by considering the 
two-point correlation function 

rW(g 2 ) = i J d*xe^ x < 0|T[v(x)4j (0)]|0 > . (3) 

The QCDSR attempts to link the hadron phenomenology with the interactions of quarks 
and gluons, which is obtained by evaluating the correlation function in two ways: an 
approximate description of the correlation function in terms of intermediate states through 
the dispersion relation, a description of the same correlation function in terms of QCD 
degrees of freedom via OPE. 

In the phenomenological side, the correlation function is calculated by inserting a com- 
plete set of intermediate states with the same quantum numbers as the tetraquark state. 
Parametrizing the coupling of the J p = 2~ tensor state to the current in term of the 
parameter fx as 

(0|j^(0)|A) = /^, (4) 

where is the relevant polarization tensor. Using Eq. (jlj) in the phenomenological side 
of Eq. (jHJ), we obtain 



n fx 



-j{g^ a 9up + 9np9ua) \ + other structures + (5) 



m x — q 2 1 2 

where the only structure which contains the contribution of the tensor meson has been 
kept. In calculations, we have performed summation over the polarization tensor using 



Efiu£ a f3 — ^T^aTyf} + ^TfipTva ^T^ u T a /3, (6) 
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with 

T^ u = -g^ + ^f 1 . (7) 

In the OPE side, we single out the structure \{g m g v p, + g^gva) whose coefficient is 
denoted as 



n «(g 2 )= [°° ds 

J Ami. 



p OPE (s) 



Ami S — q 2 



where the spectral density is p (s) = ilmn( 1 )fs). After equating the two sides, assum- 
ing quark-hadron duality, and making a Borel transformation, the sum rule can be written 



as 

J Ami 

with M 2 the Borel parameter. 

In calculations, we work at the leading order in a s and consider vacuum condensates up 
to dimension six, with the similar techniques in Refs. 26|, |27| • After tedious calculation, 
the concrete forms of spectral densities read: 

pOP^ s ) = p Pert (s) + p (m)( s) + p <^ 2 > (s) + p {ma. Gq){s) + p (?9> a (s ) > (10 ) 

with 



5 * 2 13 7T 6 Ja min a 4 J I3 min /3 4 

-f? + (i- l)r(m„s) 5 , 

O * Z 7T Jamin W J Pmin P 

„(*">(.) = ^g!L m? f - * f („ + „ - l)»(2a» + „ - 2fl» - 



3 3 * 2 13 7T 6 Jamin tt 4 Jf5 min (5 

(ggcr ■ gg) r*ma X da d/3 2 R , aR 2 n 2 flW n2 
^TT^ — mc / — / -52-(4a /3 + 6a/3 - 2a/3 - a - P)r(m c ,s) , 

(gqv ■ Gg) p*™* da f 1 ^ d(5 2 2 
oiFl — mc / — / tf^ 2 - 5 +2a/3-2/3-a)r(m c ,s) , 



with r(m c , s) = (a + (3)m 2 — af3s. The integration limits are given by a m i n = [1 
yfl - 4m 2 /s s )/2, a max = (l + yfl - 4m 2 /s s )/2, and (3 min = am 2 c /(sa - m 2 c ). 
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TABLE I: Upper limits in the Borel window obtained from the sum rule for different values of 



(GeV) 


Mi ax (GeV 2 ) 


4.6 


2.3 


4.7 


2.4 


4.8 


2.6 


4.9 


2.8 


5.0 


2.9 



To extract the mass mx, we take the derivative of Eq.® with respect to and then 
divide the result by itself 

m \ = f S ° dsp OPE {s)se- s / M2 / f° dsp OPE {s)e- s / M2 . (12) 

Before the numerical analysis of Eq.f fT2|) . we first specify the input parameters. The 
quark mass is taken as m c = 1.23 GeV 28j. The condensates are {qg) = — (0.23) 3 GeV 3 , 



(gqa ■ Gq) = m 2 (qq), m 2 = 0.8 GeV 2 , and (g 2 G 2 ) = 0.88 GeV 4 [IS]. Complying with 
the standard procedure of the QCDSR, the threshold so and Borel parameter M 2 are 
varied to find the optimal stability window. There are two criteria (pole dominance and 
convergence of the OPE) for choosing the Borel parameter M 2 and threshold Sq. In general, 
the continuum threshold Sq is a parameter of the calculation which is connected to the mass 
of the studied state, by the relation y/so ~ {m x + 0.5 GeV). 

Concretely, the contributions from the high dimension vacuum condensates in the OPE 
are shown in FigflJ We have used y/so > 4.6 GeV. From this figure it can be seen that 
for M 2 > 2.0 GeV 2 , the contribution of the dimension-6 condensate is less than 16% of 
the total contribution and the contribution of the dimension-5 condensate is less than 20% 
of the total contribution, which indicate a good Borel convergence. Therefore, we fix the 
uniform lower value of M 2 in the sum rule window as M^ in = 2.0 GeV 2 . The upper limit of 
M 2 is determined by imposing that the pole contribution should be larger than continuum 
contribution. Fig. |2] demonstrates the contributions from the pole terms with variation 
of the Borel parameter M 2 . We show in Tabled the values of M^ ax for several values of 
y/so. In FigJHJ we plot the tetraquark state mass in the relevant sum rule window, for 
different values of y/sQ. It can be seen that the mass is very stable in the Borel window 
with the corresponding threshold ^/sq. The final estimate of the J p = 2" tetraquark state 
is obtained as 



m x = (4.38 ± 0.15) GeV. 



(13) 
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FIG. 1: The OPE convergence for the J p = 2~ tetraquark state. The I and II correspond to the 
contributions from the D = 6 term and the D = 5 term, respectively. Notations a, (3, 7, A and 
p correspond to threshold parameters y/so = 4.6 GeV, 4.7 GeV, 4.8 GeV, 4.9 GeV and 5.0 GeV, 
respectively. 
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FIG. 2: Contributions from pole terms with variation of the Borel parameter M 2 in the case 
of J p = 2~ tetraquark state. Notations a, (3, 7, A and p correspond to threshold parameters 
^ = 4.6 GeV, 4.7 GeV, 4.8 GeV, 4.9 GeV and 5.0 GeV, respectively. 

In summary, by assuming X(3872) as a [eg] [eg] tetraquark state with quantum numbers 
J p = 2~, the QCDSR approach has been applied to calculate the mass of the resonance. 
Our numerical results are mx = (4.38±0.15) GeV, which indicates that X(3872) is unlikely 
to be a J p = 2~ tetraquark state. Thus, J p = 1 + assignment for the quantum numbers of 
the X(3872) is favored. 
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FIG. 3: The mass of the J p = 2~ tetraquark state as a function of M 2 . Notations a, /3, 7, A and 
p correspond to threshold parameters y/so = 4.6 GeV, 4.7 GeV, 4.8 GeV, 4.9 GeV and 5.0 GeV, 
respectively. 
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